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Abstract 

Exact solitary wave solutions of the one-dimensional quintic complex 
Ginzburg-Landau equation are obtained using a method derived from 
the Painleve test for integrability. These solutions are expressed in 
terms of hyperbolic functions, and include the pulses and fronts found 
by van Saarloos and Hohenberg. We also find previously unknown 
sources and sinks. The emphasis is put on the systematic character of 
the method which breaks away from approaches involving somewhat 
ad hoc Ansatze. 
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1 Introduction 



A number of non-integrable, non-linear, dissipative partial differential equa- 
tions (PDEs) are known to display a wide variety of complex behavior, where 
the global time evolution is often governed by the dynamics of spatially 
localized structures. For example, in defect-mediated turbulence 0], the 
disordered creation, motion and annihilation of topological defects play a 
prominent role. In some cases of spatiotemporal intermittency, well-defined 
localized objects "carry" the disorder and act as spatial delimiters of laminar 
regions 0]. This was shown in particular for the Nozaki-Bekki family of ex- 
act solutions of the one-dimensional supercritical complex Ginzburg-Landau 
equation (see below, Eq. (^) [|]. More strongly disordered regimes of simple 
PDEs like the Kuramoto-Sivashinsky equation, in which localized objects do 
not appear in an obvious manner, have been argued to be well-described 
as a "gas" of interacting pulses which are themselves exact solutions of the 
governing equation p. 

All these objects are, to a variable extent, reminiscent of the solitons 
of completely integrable equations. They can be thought of (at least on a 
qualitative level) as the dissipative counterparts of invariant tori in chaotic 
Hamiltonian systems, in the sense that they represent the preserved part 
of the rich mathematical structure of nearby integrable systems. Methods 
developed to investigate integrability of differential equations can therefore 
be expected to shed new light on spatiotemporally chaotic dynamics 



Following this line of thought, we present new exact particular solutions 
for the quintic Complex Ginzburg-Landau (CGL) equation in one spatial 
dimension, using techniques derived from the Painleve test for integrability 
[0]. The quintic CGL equation reads: 

8 A 8"^ A 

— = eA + ib,+ ic,)^ - ih - tCs)\A\'A - (b, - ^c,)\A\'A, (1) 

where e, bi, ci, 63, C3, 65, C5, are real constants and the field A{x, t) is complex. 

Eq. (^ is a one-dimensional model of the large-scale behavior of many 
nonequilibrium pattern-forming systems When the real parameters 63 
and 65 are respectively negative and positive, Eq. (|ip accounts for the slow 
modulations in space and time of an oscillatory mode close to a subcritical 
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Hopf bifurcation. The discontinuous character of this symmetry-breaking 
bifurcation is responsible for the occurrence of metastable states separated 
by fronts Pulse solutions have also been argued to exist and play an 
important dynamical role |]TD[. Examples of relevant experimental contexts 



include binary fluid convection |TI| and Taylor- Couette flow between counter- 
rotating cylinders 

We build on the work of W. van Saarloos and P.C. Hohenberg who 
recently reviewed the properties of the solutions of Eq. ([^) which they called 
coherent structures in order to emphasize their strong, usually exponential, 
spatial localization. In order to attempt to stop the proliferation of notations, 
we have chosen to use their notation throughout this paper, as well as their 
vocabulary to distinguish among the various types of coherent structures. 

Exact, analytical solutions of the CGL equation are scarce 1^, 14], and 



in any case limited to the uniformly propagating case, i.e. solutions of the 
form: 

A{x,t) = 6-'"^' A{^ = x-vt). (2) 

The original PDE, depending upon (x, t), is thus reduced to a second-order 
ordinary differential equation (ODE) in the ^ = x — vt independent variable, 
where w is a constant velocity, and " ' " stands for differentiation with respect 
to ^: 

- vA' = {e + iu)A + {hi + ici)A" - (63 - ^cg) - (65 - ic^)\A\^A. (3) 

Eq. possesses two types of fixed point: linear fixed points which cor- 
respond to a trivial vacuum state, A{x^t) = 0, and non-linear fixed points 
which correspond to plane waves: 

A{x,t) = a;ve"^^^*+''?^^ 

The next step consists in looking for connections between any two of these 
elementary objects in the phase space of Eq. (^. These connecting objects, 
the coherent structures mentioned above, are classified as follows: 

- pulses, i.e. homoclinic orbits between two vacuum states (linear fixed 
points) at ^ = ±00, 

- fronts, i.e. heteroclinic orbits between a vacuum state at ±00 and a 
plane wave (non-linear fixed point) at =Foo, 
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- sources, i.e. heteroclinic orbits between two outgoing waves at ±00, 

- sinks, i.e. heteroclinic orbits between two incoming waves at ±00. 

Sources and sinks can be distinguished from each other by checking the signs 
of the group velocities of the asymptotic plane waves at ±00. The multi- 
plicity of these coherent structures admits an upper bound, determined by 
the dimensionality of the connecting manifolds, flowing from the unstable 
eigendirection(s) of one fixed point into the stable eigendirection(s) of an- 
other fixed point. 

Exact solutions to Eq. were obtained in by using an ad hoc, re- 
duction of order Ansatz, where the first-order derivatives of the phase and 
amplitude of A{C,) are given a priori expressions. Van Saarloos and Hohen- 
berg found exact pulses and fronts, but did not mention the existence of any 
exact source or sink for the quintic equation, although such solutions are 
allowed by the counting arguments. Finally, they gave numerical evidence 
of the important role played by these special, highly non-generic solutions, 
which were shown to be "dynamically selected" in certain regions of param- 
eter space for sufficiently localized initial conditions. 



The purpose of this paper is to show that a coherent mathematical frame- 
work can foster a better understanding of these exact solutions and of the 
precise reasons for their functional form. The emphasis is laid upon a local 
study of the analytical structure of possible solutions in the complex plane, 
as opposed to the more geometrical methods referred to above. Similar ar- 
guments were recently used by Conte and Musette in the case of the cubic 
CGL equation |jl5|, for which the quintic term vanishes: 
OA d'^A 

— = eA + {h + ici)-^ - (&3 - ^C3)\A\^A. (4) 



All the known solutions of (^) were naturally retrieved in |]T5[, including the 
one-parameter family of sources originally due to Nozaki and Bekki . 

In addition to the known pulses and fronts of we find for the quintic 
CGL equation a new set of sinks and sources, whose existence is restricted 
to a low co-dimension subspace of the full (e, bi, Ci, 63, C3, 65, C5) parameter 
space. These three types of solution locally obey the same singularity struc- 
ture, and use hyperbolic functions as the elementary units from which their 
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global functional form is built. However, their respective multiplicities do 
not reach the upper bounds set in 



2 Methodology 

2.1 The Painleve test for integrability 

Our guideline will be integrability in the sense of Painleve: an ODE will be 
called integrable if its general solution is free from movable critical points. 

Let us first define these adjectives. A critical point (of a complex-valued 
application) is a point around which several determinations of the application 
occur. Examples include algebraic and logarithmic branch points. 

A movable singular point (of a solution of a DE) is a singular point 
whose location in the complex plane is not determined by the coefficients of 
the DE. This location can only be obtained from the initial conditions of the 
differential problem, i.e. from integration constants. The simplest example 
is provided by the Bernoulli equation: 



admits a movable simple pole at Xq. Conversely, singular points whose lo- 
cation depends only upon the coefficients of the DE are called fixed. Linear 
DE's only admit fixed critical points. 

In this context, integrability is intimately connected with single- valuedness: 
integrating a differential equation is ultimately equivalent to expressing its 
solution in terms of functions, i.e. single- valued applications of C into C. 
Mult i- valuedness, expressed through the occurrence of critical points, can 
be easily dispensed with when the critical points are fixed, for instance by 
removing from the domain of the solution a line in C between the critical 
point and a point at infinity. On the other hand, movable critical points are 
sources of persistent multi-valuedness, and therefore preclude integrability. 

Implicit in this scheme is the necessity of extending the domain of inde- 
pendent variables to the complex plane. Although unphysical at first sight, 




whose general solution 



u{x) 



1 
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this prerequisite is simply analogous to solving real algebraic equations for 
complex unknowns. 

The strength of the Painleve test for an ODE 



E 



(5) 



lies in its easy, algorithmic implementability. Its main requirement is the 
existence of all possible solutions u{x) of Eq. (|^) expressed as a Laurent 
expansion in a neighborhood of a movable singularity xq: 

oo 

u{x) = xi^y x{xy, (6) 

j=0 

where a is the leading-order exponent, x{^) the expansion variable, and 
{uj, j > 0} a set of constant coefficients. Following the invariant formula- 



tion of Conte |jT8|, we distinguish here the expansion variable xi^) from the 
singular manifold of Weiss, Tabor, and Carnevale 0, and only require x{^) 
to behave as a simple zero near xq: x(^) ^ x — xq. This expansion is substi- 
tuted into E[d/dx,u{x)]. The u/s are determined from recursion relations 
that develop when the coefficients at each order of x are required to vanish. 
The leading-order exponent a is determined by equating the exponents of 
the dominant order terms in the DE (|^). 

Necessary conditions for an ODE to pass the Painleve test are: 

1. the leading order a is an integer, 

2. the recursion relation for the coefficients uj can be consistently solved 
to any order. 



and possibly some other conditions not detailed here This procedure 

checks that the Laurent-type expansion for u{x) (Eq. (|^)) is both consistent 
and free from logarithmic branch points. 



2.2 Painleve analysis for nonintegrable equations 

The general solution of nonintegrable equations will fail the Painleve test at 
one of these two steps. However, this does not forbid the existence of partic- 
ular solutions, provided that they respect the singularity structure derived 
from the leading-order analysis. 
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The next step consists in determining how the local, analytical structure 
valid in a neighborhood of a singular point xq can be taken into account 
to yield global results, namely expressions of u{x) valid for C as a whole. 
Once again, the guideline is provided by one of the early results of Painleve 
on integrability the solutions of all known integrable non-linear ODEs 



of order at most two and degree one in u"{x) can be expressed as linear 
combinations of logarithmic derivatives of entire functions, whose coefficients 
are entire functions |^|. We mention here as an example the case of (P2), 



one of the six integrable second-order equations of Painleve: 

u" = 2u^ + xu + a 

where a is a constant coefficient. Its general solution can be expressed as: 

u{x) = S^LogV"! - 9^LogV'2, 

ipi and tp2 being two entire functions. 

This result should not come as a surprise, since dxljogip is by construction 
single- valued, and behaves like a pole. In this respect, the integrable (thus 
single-valued) part of nonintegrable equations is naturally expected to be 
expressible in terms of logarithmic derivatives of entire functions. 

We now turn to the definition of the class of possible solutions we consider. 
Arguments will remain mostly heuristic, although our Ansatz can be derived 
within a rigorous mathematical setting, taking into account the inherent 
invariance of Painleve analysis under the group of homographic transforma- 
tions. For more details, we refer the more mathematically-oriented reader to 



the articles |jT5[ and [|T8| and to the lecture notes pl[| . 



2.3 Ansatz for the quintic CGL equation 

Leading-order analysis for Eq. @ is achieved by balancing the highest-order 
derivative with the strongest nonlinearity. A(^) being a complex field, this 
must be done by writing two complex conjugate equations for A(^) and 
B{^) = A*{^), where "*" denotes complex conjugation. The fields A and 
B are now formally considered as independent variables, and obey: 



{bi-tCi)B" ~ {b5 + tC5)A^B\ 



2R3 (7) 



7 



Using x(0 ~ ^ ~ ^0 and x'iO ~ and feeding the leading-order Ansatz 
A Aq x"", B Bo x'^ into Eqs. (g) leads to: 



A ~ Aox~ 
B ~ Bo X' 



2 I «"U 
1 



(8) 



where Ao, Bq and cto ai'e solutions of the equations: 

J + 27^ ao - fx = 0, 
(Ao5o)2 = 2 «o, 

and the intermediate variables TZ and X are defined as: 

n = Re [ (61 +ici)(65 + 205) ] =6165-0105, 
X = Im [ (61 + ici)(65 + 205) ] =6105 + 0165. 

Without loss of generality, we assume that Aq = Bq are real constants in the 
rest of this paper. 

We consider here the non-degenerate case X 7^ 0, where: 



(9) 



The degenerate case X = 0, ao = is treated in the Appendix. It should 
be noted that all results we present (including the Appendix) respect the 
leading-order balance (|^), and are thus valid only when the coefficients of 
the quintic term and of the highest order derivative are both non-zero: 

I ^1 + ^ 0' (10) 
1 65 + ^05 ^ 0. ^'^^ 

The leading-order exponent — | + iao is not an integer: the nonintegrable 
complex Ginzburg-Landau equation has already failed the Painleve test. As 
expected, its general solution cannot be expressed in terms of elementary 
functions. However, partial integrability remains possible, in so far as one 
looks for solutions exhibiting minimal - but necessary - multi-valuedness, 
including movable algebraic and logarithmic branch points. This requirement 
is clearly fulfilled by the following rewriting of Eq. (§) : 

A(0 = Ao e"*^* e^°°®(«). (11) 
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We introduced here a generalized amplitude -R(0 and a generalized phase 
0(^), a priori complex- valued. The amplitude R{C) is assumed to have at 
worst the singular behavior of a pole, that of a logarithmic branch point. 

The global structure of possible solutions is introduced as follows. To a 
given leading-order exponent correspond four distinct values of Aq, related 
through a phase shift of |. Four families of entire functions ipi, i = 1, ...,4 
must therefore be introduced in the expression of R{C)- 

R{0=ro{0+^2^^i0^^^ogiJi. (12) 
1=1 

The coefficients r^, i = 0, ...,4 are assumed to be entire functions of ^, thus 
ensuring a pole behavior for R{C,) in the vicinity of any of the movable zeroes 
of the ipi^s. For simplicity, we will restrict our Ansatz to the case where: 

- only two families ipi and ip2 are used, and 

- the coefficients tq, ri and r2 are real constants, 

thus keeping the number of unknowns at a tractable level with respect to 
the number of equations. Leading-order analysis leads to ri = ±r2 = ±1, 
when conducted according to Eq. (|T^) with ro, ri, r2 real constants, and 
rs = r4 = 0. Up to a constant phase shift, we can fix ri = +1, r2 = ±1. Our 
Ansatz for R{C,) reads: 

RiO = ro + S^LogV'i ± a^LogV'a. (13) 

Let us now turn to the Ansatz for 0(.^). Only one family is necessary here, 
which we denote ipi. Since the complex Ginzburg-Landau equation is invari- 
ant under an homogeneous phase translation A A e'^'^, only contributes 
through its gradient d^Q. We only need to define an expression for d^Q, with, 
again, a pole singular behavior at worst. This expression reads: 

a^e(o = ^o + 5?Log^i, (14) 

where a constant coefficient 9o was introduced. 

The entire functions ipi and ip2 are next defined as solutions of integrable 
differential equations. For simplicity, we assume that the ■z/'j's are solutions 
to the second-order linear ODE: 



de 4 



(15) 
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Linear first-order differential equations have been excluded. Their solutions 
have constant logarithmic derivatives, leading to trivial expressions of -R(C) 
and 0(0 which correspond to the fixed-point solutions of Eqs. (|^). Other 
choices of the defining ODE are possible in principle, but were at first dis- 
carded, again for reasons of simplicity. This helped in keeping the algebraic 
manipulations needed later at a tractable level. 
The general solution to Eq. (|15|) reads: 

V'(0=V^o cosh^(e-eo), (16) 

where ipo and are the two integration constants. The value of ipo can be 
set to 1, since only the ratio d^ipi/ilji contributes. We define the two families 
as two independent solutions to Eq. (p!6[) separated by a constant phase shift 
denoted ka: 

= cosh 1(^-^0 + a) 



^2(0 = cosh 1(^-^0 -a 
Elementary manipulations then show that [p!5|: 



(17) 



cosh/c(4— 4o) + cosh[ka) 
9,Log^, - ^^Log^2 = eoshfc(^-eo) + cosh(fca) - ^ 



Our Ansatz (|T^ can be expected to yield all possible solutions to the com- 
plex Ginzburg-Landau equation involving hyperbolic functions or, up to lin- 
ear combinations, exponentials. Hyperbolic tangents can be obtained from 
Eq. (|18D when cosh(fca) = 0, i.e. ka = i|, hyperbolic secants from Eq. (|1^) 
when cosh(fca) = 0. 

Eq. (|l^) shows that two families having identical coefficients are in prac- 
tice equivalent to one family, up to dividing the wavenumber /c by 2 and to 
setting ka = i^. This allows a slight modification of the equation defining 
0(0, formally using the two families ipi and ijj2. The complete Ansatz for 
A{^) now reads: 

A(0 = Ao e"*"* i?(0^ e^°o®(«), 

R = ro + S^LogV^i ± 9^Log'02, /20) 
d^Q = 9o + di^Logijji + 9^LogV'2, 

Albeit very restrictive, this Ansatz suffices to retrieve all the exact solutions 
quoted in Jll. 
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2.4 Computational aspects 

From now on, the constant coefficients 6q, uj, c will be real numbers, in order 
to avoid unbounded solutions. The wavenumber k is also assumed to be real. 
Possible periodic solutions are thus discarded, since they lack the required 
spatial localization. On the other hand, taking ro G C would not restrict 
the generality of Ansatz (|20|) . However,we checked that the quintic equation, 
unlike the cubic one (see [0), does not admit any solutions respecting the 
less restrictive hypothesis tq = xq + i yo; Xq, z/q G R; XoUq ^ 0. 

The basic principle of our resolution is to turn the differential equation 
we want to solve into a much simpler, purely algebraic problem. This is 
made possible by the analytic considerations of the previous section: the 
spatial structure of solutions to Eq. (|^) is supposed to be fully contained 
in the elementary logarithmic derivatives S^Log^/'j, whose functional form 
will never be made explicit in our computations. The sometimes intricate 
algebraic manipulations can then be solved quite easily with the help of any 
symbolic mathematics package, such as Mathematica |^, or AMP p3[ . 

We ffist substitute our general Ansatz (Eqs. (|20|) ) into Eq. (H) and suc- 
cessively eliminate all derivatives of ^/'j of order greater than or equal to two 
by using: 

de^ = T 

Eq. d^) is then equivalent to a polynomial equation in the S^LogV'i and 
9^Log'02 variables: 

E E (5«Log^i)"^ (5cLogV^2)" = 0, (21) 

fe=0 m+n=fc 

where the coefficients depend algebraically on the parameters (e, 6j, Cj), 
on the unknowns tu. A;, f , a, tq, 6^0, and on and Oq, whose values are known 
from Eq. (H). 

A convenient way of taking the phase shift ka into account is to use a new 
variable /iq, defined as /io = coth(fca). The products (9^Log'i/'i)"^ ((?^Log?/'2)" 
can be recursively linearized, from m + ?T, = 4tom + ?T, = 2, by means of the 
following identity: 

k^ k 

d(Logi)x d^Logilj2 = Y ~ /^o 2 [ ^^Log^i - 9^LogV'2 ] • (22) 
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The coefficient /io is now a real constant, in agreement with the previous 
assumptions. We obtain: 

j^E-j {d^Log^iY + Eo + j^Ej {diLog^2y =0. (23) 

Solving Eq. ( [23| ) amounts to canceling all the Ej coefficients, j = —4, . . . , 4. 
This is done recursively, from E_4 and £^4 to E'o, by a triangularization tech- 
nique. Parameters and unknowns are considered on an equal footing, as 
variables whose values are successively obtained from Ej and substituted 
into equations Ei, < |i| < |j|, thus decreasing the number of unknowns in 
equations of lower index. For each j, the pivoting variable was determined 
by singling out which variable admits the simplest expression with respect to 
all other variables, while excluding possibly vanishing denominators. These 
selected variables are: 

- as = 65 — ics, obtained from i?_4 or E4, 

- = bs — ic3, obtained from -E_3, 

- i) = V + 2iaoai9o, obtained from E3, where ai = 61 + ici, and 

- e = e + iOouj, obtained from E^2- 

Equation E2 is then seen to vanish, and the remaining E_i, Eq and Ei 
equations depend only on the parameters bi, Ci, the leading-order quantities 
Aq, ao, and the unknowns k, /xq and vq. Systematic resolution of these 
equations lead to the three cases detailed below, where the explicit values 
of the real unknowns k, v, uj, 6q, vq and /iq are given as functions of the 
parameters {e, bi, Ci), Aq and ao- The leading-order quantities Aq and ao are 
considered as parameters, since their value can be expressed as functions of 
61, Ci, 65, and C5 (Eq. (D). 

3 Results 

We do not mention here "unphysical" solutions obtained for complex values 
of the parameters k, v, u, 9q, tq or /iq, although such solutions may well 
become "physical" when more general hypotheses are considered. These 
questions have been left for future work. 
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3.1 Pulses 

The pulse solutions given in can be obtained from the following Ansatz: 

i R = S^LogV^i - d^Log'ip2, 

I = 6*0 + S^LogV'i + d^Logip2, 

where ro was set to in order to respect the suitable asymptotic behavior: 



The solution reads: 



A{x,t) 



lim A{^) = lim A{^) = 0. 



j^^^-iu^t^iaoeo^ [ cosh - ^o) + cosh ka f 
k sinh(A;a) 



(24) 



cosh — (^o) + cosh(A;a) 



In this particular case, triangularizing the equations Ej, j = —4, ...,4 leads 
to -0 = 0, or t> = —2iao{bi + ici)9o. The velocity v being a real constant, we 
obtain: 



bi0o = 0, 

V = 2ciao9o. 

We find a set of solutions, discrete when bi ^ 0, v = 9q = 0, parametrized 
by the velocity v or by Oq when 6i = 0, in a co-dimension-one subspace of 
parameter-space defined by: 



C3 



— 2al) + 3aoCi 



Saobi + ci{2al - 1] 



The unknowns k, fiQ and uj can be computed from: 

4£ 



k' 



kfio 



ijj 



- 4ag) +4aoCi' 



&i(l -2^2 

-CiOo^'o 



o; -r SttoCi 

e — 4ao^i + Ci(l — 4ao) 
ao — 4ao) + 4aoCi 
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3.2 Fronts 

The front solutions of involve hyperbolic tangents. The Ansatz we use 
goes as follows: 

( R = To + ^^LogV'i + ^^Log^s, 

[ d^Q = 60 + SgLogV'i + 9gLogV'2- 

Solutions were found only when /iq = 0. Their explicit form reads: 

A{x,t) = Aoe-*^*e^°°^°« [cosh A;(e - ^o)]*"" [A;(tanhA;(e - ^o) ± (25) 
where the appropriate asymptotic behavior is obtained by setting ro to ±k: 
hm A{0 = 0, 

lim A(0 = y±2fcAoe-^'^*e^"«(^°^'=\ 

Its variables obey different relationships according to the value of bi: 

- First subcase: 61 7^ 0. A discrete set of solutions is found. The 
wavenumber k is determined by a quadratic equation: 

Ofcfc^ + bkk + Cfc = 0, 

where: 

-(1 + 4^2)2 (36? + 4c?), 

T2AI (1 + Aal) {bibs - 2ciC3 + 2ao(6iC3 + 263C1)), 
bie (1 + 4a2)2 - (4 + ib^al - 4006303). 

The unknowns 6*0, v and lj are obtained as functions of k: 
^ - ±kil-2^) + Al ^3-2«o&3 



'aobi ao(l + 4ao)^i' 

r, bl + 4 u o ^0 (C1C3 - 6163) + 2ao(6iC3 + 63C1) 
bi bi 1 + Aal 



6f«o(l + 4«^)2 



where: 



= -Ci(l+4aoT(56? + 4c?), 
b^ = ±2Al{l + 4al) 

X [{bjcs + 2cic3 - 26163C1) - 2ao {b% + 263c? + 261C1C3) 

= (C3 - 2^063) [(26163 - C1C3) + 2^0(6301 + 261C3)]. 
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Second subcase: bi — 0. In the non-degenerate case treated here 
{I — 61C5 + C165 7^ 0), this condition imphes that ci cannot vanish. We 
find a one-parameter family of solutions restricted to a co-dimension- 
two subspace of parameter space defined by: 

bi = 0, 

ci{l + Aalfe = {c3-2aob3){b3 + 2aoC3) A^. 

The unknowns ^0, and k are parametrized by the velocity v: 

, Al 263(1 -a2) + 5c,oC3 



^0 = 



UJ — 



V 

2q;oCi 2aoCi 



AaoCi 



At 



4q;o(1 + 4:aiyci 



1 + 4ag 



[(263 + 4aoC3)'- (2ao?^3-C3)2], 



Al C3-2ao&3 



2ci l + 4:al 
3.3 Sources and sinks 

Sources and sinks are obtained from the Ansatz: 

' R = To + S^LogV'i - d^Log'ijj2, 
< d^Q = 9o + d^Logipi + d^Logip2, 
.ro 7^ 0. 

Solutions exist for /Iq ^ only, and interpolate between two plane waves: 

lim A(0 = Aoroe-*'^*e^[°^°(^°-'=)^l, 
^— »— 00 

lim A{C) = Aoroe-''^Vl"°(^°+'=)^l. 
Their expression reads: 

A{x,t) = Aoe-^'^V'^°^°«[coshA;(e-6)+coshA;a]^"° (26) 

_ 1 

k sinh(A;a) 



' cosh A;(^ — ^0) + cosh(/ca) 
As for pulses, we obtain v — 0, whence: 

61^0 = 0, 

V — 2ciao9o- 



+ Tq 
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We find a set of solutions, discrete when 61 ^ 0, v — 60 — 0, parametrized 
by the velocity v when bi = 0, in a co-dimension-one subspace of parameter- 
space defined by: 

1 r- 



where: 

To = 



^=4 



^0 



X 

^0 



7ao(^'i^'3 - C1C3) + {2al - 3)(63Ci + 61C3) 
3q;o(^'i^3 - C1C3) + {2al - l)(63Ci + 61C3) , 

3q;o(6i63 - C1C3) (1 - 2q;o)(63Ci 61C3) 



The remaining unknowns are given by 

1 



a; 



1^0 



The asymptotic group velocities Vg^ ± in a co-moving frame at ^ — > ±00 are 
given by: 



where the respective asymptotic pulsation and wavenumbers are fl — cu and 
K± — ao {9q ± k). We obtain: 



1 ( duj ^duj 



2ao \ d9o dk 
V _|_ ^1 ~ 3q;o6i 
2ao 2ao 

In the generic case 61 7^ 0, v — Oq — 0, we find 

Ci — 3q;o6i 



k. 



-Vn 



2ao 



whose sign can be either positive or negative. These solutions can be either 
sources or sinks. 
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4 Conclusion 



The search for the functional form of the coherent structures appearing in 
disordered regimes of extended, non-equihbrium systems can be made sys- 
tematic by following two basic principles: 

- the singularity structure of possible solutions of the relevant DE's must 
be taken into account at an early stage, 

- logarithmic derivatives of entire functions are the elementary units from 
which exact solutions can be built. 

In this work, we followed these principles and derived systematically an 
Ansatz for coherent structure solutions of the one- dimensional quintic com- 
plex Ginzburg-Landau equation. In addition to the already known pulses 
(Eq. (0)) and fronts (Eq. (pS])), we found new source and sink solutions 
(Eq. (0)). The corresponding computations involved only the simplest pos- 
sible Ansatz compatible with our framework. We reduced the number of fam- 
ilies of entire functions to two, used real coefficients throughout, and chose 
to define the elementary entire functions from the simplest available ODE. 
Within these restrictions, exponentials appeared as the simplest building 
blocks from which exact solutions can be formed. Relaxing these constraints 
can be expected to yield new - and more complex - solutions, provided that 
the resulting algebraic computations are not made intractable by increas- 
ing the degree of equations to be solved. These less restrictive Ansatze are 
currently under investigation. 

The stability and dynamical relevance of the solutions, especially the 



new source/sink, remain to be investigated. In |T^, a survey of the special 
analytical, topological and dynamical properties of the highly non-generic 
fronts (Eq. (]25|)) was given, and the question of relating these three aspects 
was raised. These properties are crucial steps in trying to understand the 
spatiotemporally disordered regimes exhibited by the equation, and deserve 
closer scrutiny. Another point of interest is the question of a possible role 
played by these solutions in regions of parameter space out of their domain 
of existence. There, objects related to these solutions, but out of reach of the 
simple Ansatz used in this work, may appear as the relevant building blocks 
in the (chaotic) dynamics. 



17 



To conclude, let us stress again the systematic character of the approach 
taken and the potential extensions of the case treated here, not only to less 
restrictive Ansatze, but also to other non-linear PDEs of physical interest. 



Appendix 

We treat here the degenerate case: X = 61C5 + 65C1 = 0. Leading-order 
analysis then leads to: 



«o = 0, (27) 

4 _ 3 — C1C5 
° - i hl + d ' 

whenever the condition (0) applies. This singular behavior is taken into 
account by writing: 

where -R(0 has at most the singular behavior of a pole, and 0(0 is a regular 
function. In the spirit of Section p.3| , we write: 

= To + d^{Logipi) ± a^(LogV'2), 
= Oo, (28) 

= f 2 = 1,2. 

As before, looking for particular solutions of the quintic CGL equation with 
the restriction 61C5 + 65C1 = amounts to solving a system of algebraic 
equations. Its systematic resolution leads to pulse, front, source and sink 




solutions. The method is identical to that presented in Section p^ , and 
solutions are given within the restrictions of Section ^. 

Pulses 

The pulse solutions read: 

k sinh(fca) 



cosh/c(^ — ^0) + cosh(A;a) 
and necessarily respect ■0 = 0, bi9o = 0. We distinguish the two cases: 



(29) 
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- First subcase: 61 7^ 0. A discrete set of stationary pulses is found in 
a co-dimension-two subspace of parameter space defined by: 

61C5 + 65C1 = 0, 
&1C3 + &3C1 = 0. 

Such restrictions include the case of the Real Ginzbug-Landau equation 
(RGL: ci = C3 = C5 = 0): 

hlAl^A-b^lAl^A. (30) 

0, 
0, 

-As, 

-hAl 

ec\. 

- Second subcase: bi = 0. A two-parameter family of pulses is found 
for the Quintic-Cubic Schrodinger equation (ci 7^ 0, C5 7^ 0): 

The free parameters are chosen to be the wavenumber k and velocity 
v: 

2ci6o = V, 

Aciuj = —v"^ — cfk'^, 

cikfiQ = A^cs. 

Fronts 

Fronts are obtained for Ho = and tq = ±k: 

= A;^>loe^'^°^"'"*^[tanhA;(e-6) ±1]^- (32) 
The condition 61 + ici ^ leads us to distinguish two cases: 

- First subcase: 61 7^ 0. A discrete set of fronts is found in a co- 
dimension-one subspace of parameter space defined by: 

61C5 + 65C1 = 0. 



eA + b ^ 
dt ^ dx^ 



All parameters are fixed: 



V 

Oo 

bik^ 

bik/io 
biCU 
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The wavenumber A; is a solution of the quadratic equation: 

dkk'^ + hkk + Cfc = 0, 

where: 



= (36? + 4c?), 
= ±2Al (6163-2C1C3), 
= -ehi + clA^. 

The unknowns 6*0, v and oj are obtained as functions of k: 



biOo = T'2cik +Alc3, 
biv = T2{bl + cl)k + 2Al{ciC3-hh), 
h\uo = a^P + b^k + 



where: 



= -ci(5fe2 + 4c?), 
= ±2A2(62c3 + 2c2c3-26i63Ci), 
= C3 (26163 - C1C3). 

- Second subcase: bi = 0, ci 7^ 0. A one-parameter family of solutions 
is found in a co-dimension-three subspace defined by: 

bi = 0, 
&5 = 0, 
eci = A'^bscs, 

thus including the generalized quintic-cubic NLS equation (Eq. (|3TD) 
when e = 63 = 0. All coefficients can be expressed as functions of the 
velocity v: 

2ci^o = v + 2Albs, 

icico = -v^ + A^{4bl- C3), 

2cik = ±Alc3. 

Sources and sinks 

These solutions read: 

A(x,t) =Aoe't^o^-"*l ro + 



k sinh(A;a) 



cosh /c(^ — ,^0) + cosh(/ca) 
and respect -0 = 0. We distinguish the two cases: 



(33) 



First subcase: 61 7^ 0. A discrete set of stationary solutions is found 
in a co-dimension-two subspace of parameter space defined by: 

61C5 + 65C1 = 0, 
61C3 + 63C1 = 0, 
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including the RGL equation (Eq. (^)). All parameters are fixed: tq is 
determined by the quadratic equation: 



3blrl + AAlb^ro -Ae = 0, 

and: 

V =0, 
^0 = 0, 

6iP = -2Alhro + 4e, 
biu = eci, 

The far-field group velocity Vg vanishes, due to an asymptotic stationary 
wave behavior: 



lim A{x,t) = Aoy/r^e' 



-iujt 



Second subcase: 6i = 0. A two-parameter family of sources or sinks 
is found for the quintic-cubic NLS equation (0). We choose to use v 
and Vq as free parameters (ci 7^ 0): 

2ci9o = V, 

Aciu = -v^ + 3c^rg - 4AgciC3ro, 
= 3cir2 - 2Agc3ro, 

/^O ~ 2fcro ■ 

The asymptotic group velocities Vg^ ± read: 
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